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Abstract 

The Tasmanian Devil Facial Tumour Disease (DFTD) provides a unique oppor- 
tunity to elucidate the long-term effects of natural and anthropogenic selection 
on cancer evolution. Since first observed in 1996, this transmissible cancer has 
caused local population declines by >90%. So far, four chromosomal DFTD vari- 
ants (strains) have been described and karyotypic analyses of 253 tumours 
showed higher levels of tetraploidy in the oldest strain. We propose that increased 
ploidy in the oldest strain may have evolved in response to effects of genomic 
decay observed in asexually reproducing organisms. In this study, we focus on 
the evolutionary response of DFTD to a disease suppression trial. Tumours 
collected from devils subjected to the removal programme showed accelerated 
temporal evolution of tetraploidy compared with tumours from other popula- 
tions where no increase in tetraploid tumours were observed. As ploidy signifi- 
cantly reduces tumour growth rate, we suggest that the disease suppression trial 
resulted in selection favouring slower growing tumours mediated by an increased 
level of tetraploidy. Our study reveals that DFTD has the capacity to rapidly 
respond to novel selective regimes and that disease eradication may result in 
novel tumour adaptations, which may further imperil the long-term survival of 
the world's largest carnivorous marsupial. 



Introduction 

More than 35 years ago, NoweU (1976) suggested that can- 
cer progression should be regarded as an evolutionary 
process. We now know that cancer is subjected to selective 
regimes similar to those experienced by asexually reproduc- 
ing organisms (Merlo et al. 2006). Cancer cells, like other 
asexual organisms, do not undergo meiotic recombination. 
How tumour cells survive the loss of heterozygosity and 
the emergence of recessive mutations remains unresolved 
(Sole and Deisboeck 2004). Proposed mechanisms to coun- 
teract genomic decay include chromosomal rearrangements 
that alter normal cell cycles and apoptotic responses, chro- 
mosome breaks and tolerance of deleterious mutations 
(Merlo et al. 2010). 



The presence of aberrant karyotypes in malignant cells 
was first observed over a century ago (von Hansemann 
1890) and led to the recognition of the role of missegregat- 
ing chromosomes in tumour development (Boveri 2008; 
Chow and Poon 2010). Indeed, most malignant tumours 
have been found to harbour structurally and numerically 
rearranged chromosomes and multiple centrosomes, the 
conjoint causes and consequences of abnormal mitosis 
(Sen 2000; Storchova and Pellman 2004; King 2008; Storch- 
ova and Kuffer 2008; Ganem et al. 2009; Little 2010; 
Mosieniak and Sikora 2010). Apart from segmental chro- 
mosome defects and single chromosome losses cancer cells 
exhibit altered ploidy, with chromosome numbers ranging 
from hypodiploid (i.e. having a chromosome number 
lower than the diploid number) to hypertetraploid 
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(i.e. having a chromosome number greater than but not an 
exact multiple of the normal diploid number; for reviews 
see: Storchova and Pellman 2004; Storchova and Kuffer 
2008; Davoli and de Lange 2011). Aneuploidy and poly- 
ploidy have, however, been shown to provide cancer cells 
with adaptive potentials (Yuen and Desai 2008). For exam- 
ple, polyploidization may provide adaptive advantage to 
cancer cells by masking deleterious mutations (chromo- 
some losses, gene deletions and inactivating mutations) 
and ameliorating the effects of deleterious recessive alleles 
(Otto and Whitton 2000; Otto 2007; Davoli and de Lange 
2011). Moreover, tetraploidy allows tumour cells to sustain 
a higher mutation rate and may stimulate additional gen- 
ome structure modifications facilitating adaptive changes. 
The immediate effect of polyploidization is a general rise in 
cell volume, and slower development due to increased gen- 
ome size (Cavalier- Smith 1978; Gregory 2001; Otto 2007). 
However, changes in ploidy also upset the geometric 
machinery used to segregate chromosomes resulting in 
unstable genomes, rapid genomic repatterning and 
increased genetic diversity (Wendel 2000; Otto 2007). The 
increased genetic polymorphism associated with tetra- 
ploidy may promote the survival of certain polyploid cells, 
stabilize their genomic configuration and therefore fuel the 
evolution of polyploid cell populations. Thereby, poly- 
ploidy not only appears to promote tumorigenesis, but also 
steers cancer cell progression through a fitness landscape 
during cancer evolution. 

Although the evolution of neoplasm in human cancers 
occur on a timescale of years, anthropogenic selection, 
administered as medical treatments, has been shown to 
accelerate the development of novel and aggressive as well 
as drug resistant cancer strains (Merlo et al. 2010). A sig- 
nificant problem when investigating how such therapies 
may affect human cancer evolutionary trajectories is often 
the short lifespan of both tumours and patients. The Tas- 
manian Devil Facial Tumour Disease (DFTD) provides a 
unique opportunity to elucidate the long-term effects of 
natural and anthropogenic selection on cancer evolution. 
This contagious cancer was observed first in 1996 and is 
transmitted between Tasmanian devils {Sarcophilus harrisii) 
by biting during social interactions (Pearse and Swift 2006; 
McCallum et al. 2007; McCallum 2008; Murchison 2008). 
The disease generally results in death of infected animals 
within 6 months and has led to local extinctions of more 
than 90%, questioning the long-term survival of this iconic 
animal (Jones et al. 2007; McCallum 2008). Cytogenetic 
analyses have revealed that DFTD is caused by a rogue cell 
line (Pearse and Swift 2006), which originated from Schw- 
ann cells of the peripheral nerve sheath (Loh et al. 2006; 
Murchison et al. 2010). devil facial tumour (DFT) cells 
possess a highly rearranged genome, characterized by 
tumour-specific complex translocations and chromosomal 
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rearrangements (Pearse and Swift 2006; Deakin et al. 2012; 
Pearse et al. 2012). The clonal nature of DFTs has been 
supported by both large-scale genomic (Miller et al. 2011; 
Murchison et al. 2012), immunohistological (Loh et al. 
2006) and genetic analyses (Siddle et al. 2007, 2010; Belov 
2011). Although four chromosomal variants (strains) have 
been observed, exhibiting minor cytogenetic differences, 
the genome of DFT cells appears to remarkably be stable 
(Deakin et al. 2012; Pearse et al. 2012). 

In an attempt to reduce the prevalence of DFTD, infected 
devils, approximately 33%, were removed from one site, 
the Forestier Peninsula, in Tasmania between 2006 and 
2010 (Lachish et al 2010; Beeton and McCallum 2011). 
The disease eradication trial provides a unique opportunity 
to elucidate the long-term effects of anthropogenic selec- 
tion on DFTD evolution. In this study, we use karyotypic 
analysis to investigate overall temporal changes of tumour 
ploidy. Furthermore, we investigate the possible effects of 
the removal programme on cancer evolutionary trajectories 
on the Forestier Peninsula compared with other areas of 
Tasmania not subjected to anthropogenic selection. 

Methods 

Samples 

Tumour tissue samples used in the study were collected 
between 2006 and 2011 at 11 sites within the DFTD- 
affected areas of Tasmania (Fig. 1 which also provides data 
on number of samples collected at each of the 11 sites). A 
total of 253 diseased devils were analysed. Due to the trap- 
ping regimes employed, we were unable to obtain compa- 
rable number of samples from the 11 sites, making robust 
among-site comparisons unattainable. However, as the 
devil population on the Forestier Peninsula (Fig. 1) has 
been subjected to a disease suppression trial, that is, the 




Figure 1 Map of Tasmania showing the location of the 1 1 sites sam- 
pled and number of samples collected at each location. 
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removal of infected devils, we decided to investigate how 
such artificial selection may have affected tumour evolu- 
tion. Of the 149 tumours collected at Forestier Peninsula, 
148 were classified as strain 3 and one as strain 2. In con- 
trast to most of the other sites, no samples were collected at 
Forestier Peninsula in 2011. 

Cell culture 

Detailed description of DFTD cell culture and cytogenetic 
analysis has been previously described by Pearse et al. 
(2012). Briefly tissue biopsies and fine needle aspirates were 
transferred to sterile Petri dish and washed three times with 
Dulbecco's phosphate-buffered saline (Invitrogen, Carls- 
bad, CA, USA), with 0.1 mL penicillin/streptomycin solu- 
tion (Sigma-Aldrich, St Louis, MO, USA) and 1 mg/mL 
amphotericin B (Sigma-Aldrich). Solid tissues were then 
disaggregated in 3 mL of prewarmed AmnioMax C-100 
medium (Invitrogen/Life Technologies, Carlsbad, CA, 
USA) and homogenized in 3-mL syringe with an 18G nee- 
dle untU a milky single cell suspension was formed. Cancer 
cells were aliquoted into 20-mL culture flasks containing 
8 mL of AmnioMax C-100 media (Invitrogen/Life Tech- 
nologies), 0.1 mL of penicillin-streptomycin solution 
(Sigma-Aldrich), and 1 mg/mL of amphotericin B (Sigma- 
Aldrich). Cultures were incubated at 35°C. Tumour cells 
were harvested after 24-48 h in culture. 

Cytogenetic analysis 

Three hours prior to harvesting 0.1 mL of demecolcine at 
10 mg/mL (Sigma-Aldrich) was added to each culture. The 
cells were centrifuged for 10 min at 100 g. After the super- 
natant was discarded and the cell pellet was slowly resus- 
pended in 7 mL of hypotonic 0.075 m KCl and placed in a 
water bath at 37°C for 10 min; 2 mL of chilled Carnoy's 
fixative (3:1 ratio of methanol and acetic acid) was added, 
and the tubes were centrifuged for 10 min at 100 g. The 
pellet was gently resuspended in fixative and stored at 
— 2°C overnight. The following day, the cells were washed 
4x in fresh fixative and resuspended. Chromosome spreads 
were achieved by adding a droplet of the suspension onto a 
frozen microscope slide. Slides were subsequently air-dried 
and incubated at 57°C for 3 days. G-banding was con- 
ducted by treating slides with a 0.15% solution of trypsin 
(Sigma-Aldrich) for 30 s, then staining with Leishman's 
stain for 2.5 min followed by mounting with Leica mount- 
ing medium (Leica Microsystems, North Ryde, NSW, Aus- 
tralia) for analysis. G-banding analysis was performed 
using a Leica DM 2000 microscope (Leica Microsystems) 
and photographed with a Leica DFC 420 C camera (Leica 
Microsystems). Karyotypes were made, originally by hand 
and later (from 2008) using Video testeKaryo 3.1 software 
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(VideoTest, Saint Petersburg, Russia). At least 20 metapha- 
ses were analysed for each individual, and approximately, 
200 cells were examined. 

Statistical analyses 

Data were examined for normality before analysis, and 
when normality could not be achieved, nonparametric sta- 
tistics were employed. Logistic regression was used when 
investigating temporal variation in tetraploidy among the 
253 devil tumours where presence, that is, when all of 200 
cells in the metaphases were recorded as being tetraploid, 
was entered as '1' and total absence of tetraploids was 
entered as '0'. JMP, version 5.1 (SAS Institute Inc., Gary, 
NC, USA), were used in all statistical analyses. 

Results 

The large number of karyotypes used in this study made it 
possible to conduct robust tests to elucidate a possible vari- 
ation in tetraploid tumours among the four strains. Our 
results revealed a significant variation of tetraploid tumours 
when employing the total number of samples collected at 
the 11 sites as well as when omitting the samples collected 
at the Forestier peninsula (x^ = 11.7, P = 0.008, df = 3 and 

= 9.6, P = 0.02, df = 3, respectively. Fig. 2). Both analy- 
ses showed that the proportion of tetraploid tumours were 
highest in the strain 1, the oldest of the four strains (Fig. 2). 

We also observed a significant temporal increase in tetra- 
ploid tumours collected at the 11 sites from 2006 to 2011 
(logistic regression with tetraploidy as dependent and year 
as independent variable: y^ = 16.4, P < 0.0001, df = 1; 
Fig. 3A). However, when conducting the same analysis, but 
excluding the tumours collected from the Forestier Penin- 
sula, and hence restricting the analysis to the 10 remaining 




12 3 4 
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Figure 2 Proportion of tetraploid tumours recorded in the four DFTD 
strains. The numbers above the bars depict samples sizes. The two bars 
at strain 3 depict the proportion of tetraploid tumours when including 
(n = 154) and excluding (n = 7) the samples collected at Forestier Pen- 
insula. 
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sites, no temporal increase in tetraploid tumours was (A) 1 -00 - 
observed {yj = 0.15, P = 0.70, df = 1; Fig. 3B). A third 
logistic regression analysis, restricting the data to samples 
collected at Forestier Peninsula, revealed a significant tem- 
poral increase in tetraploid tumours among the Forestier 
devils (/^ = 35.0, P < 0.0001, df = 1; Fig. 3C), clearly 
demonstrating that the overall temporal increase in tetra- 
ploid tumours was caused by the increase in tetraploids at 
the Forestier Peninsula. 



Discussion 

Devil Facial Tumour Disease is a horizontally transferred 
asexuaUy reproducing clonal cell line, which during the last 
16 years have been exposed to the negative effects associ- 
ated with MuUer's ratchet, resulting in mutational melt- 
down and ultimately extinction. However, this obligate 
parasite has been able to survive and counteract the effect 
of deleterious mutations, genomic instability as well as 
being able to infect >100 000 devils (McCallum 2008). 
DFTD hence provides a unique opportunity to study can- 
cer evolution in vivo. 

Polyploidization is more common in asexual compared 
to sexual organisms (Otto and Whitton 2000) and provides 
an adaptive advantage to asexual organisms, such as cancer 
cells, by masking deleterious mutations and ameliorating 
the genomic decay process (Haldane 1933; Orr and Otto 
1994; Orr 1995; Otto and Whitton 2000; Otto 2007; Davoh 
and de Lange 2011). The higher number of tetraploids 
recorded in the oldest DFT strain (strain 1) hence provides 
a possible additional mechanism by which this asexuaUy 
reproducing obligate parasite has been able to avoid muta- 
tional decay. Increased ploidy has also been associated with 
slower tumour development and DFT cell growth rate 
(Pearse et al. 2012). During the disease suppression trial, at 
the appearance of first lesions, infected devils were removed 
from Forestier Peninsula (Lachish et al. 2010; Beeton and 
McCallum 2011). Such a selective regime could have 
favoured slower growing DFT cells ultimately resulting 
in the increased level tetraploid tumours observed at this 
location. 

In a recent study, Murchison et al. (2012) suggested that 
the unique mitochondrial DFTD lineage present on the 
Forestier Peninsula had most likely emerged due to a selec- 
tive sweep. Our results provide further evidence that the 
observed genetic and chromosomal changes at this site 
were most likely caused by selective sweep initiated by 
increased selection administrated via anthropogenic inter- 
action (ongoing removal of DFTD infected devils). Unfor- 
tunately, the selective culling of infected devils neither 
slowed disease progression nor reduced population-level 
impacts of DFTD and was therefore abandoned in 2010 
(Lachish et al. 2010). 
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Figure 3 Logistic regression analyses of temporal variation in tetraploid 
tumours. Figure (A) Depicts the variation from 2006 to 201 1 using our 
complete data set, that is, the analysis is based on all the 1 1 popula- 
tions. Figure (B) Depicts the temporal variation in tetraploid tumours 
excluding the samples collected at Forestier Peninsula. Figure (C) 
Depicts the temporal variation in tetraploid tumours collected at the 
Forestier Peninsula. 



©2013 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd 7 (2014) 260-265 



263 



Temporal variation of polyploidy in DFTD 

Transition from whole-genome duplication via aneu- 
ploidy to malignancy is a common feature of several 
human cancers (for review see Davoli and de Lange 2011). 
In spite of reducing cell proliferation rates, tetraploidiza- 
tion has also been linked to the metastatic, aggressive as 
well as drug resistant stages of certain human malignancies 
(Castedo et al. 2006; Davoli and de Lange 2011). As sug- 
gested by Davoli and de Lange (2011), such significant 
changes in tumour phenotypes is most likely caused by an 
enhanced ability of tetraploid tumours to sustain a higher 
incidence of mutations, thereby increasing the probability 
of adaptive changes and increasing the probability that 
evolving tumours will accumulate mutations needed to 
progress to a malignant state. If the increased level of tetra- 
ploid tumours at Forestier Peninsula results in the evolu- 
tion of a more malignant strain of DFTD, this may further 
imperil the long-term survival of the world's largest carniv- 
orous marsupial. Although future studies are needed to 
elucidate the connection between malignancies and tetra- 
ploidy in DFTD, our study clearly demonstrates that DFTD 
tumours are able to rapidly respond to increased selection 
and adapt to a novel selective regime. Due to the observed 
low genomic (Miller et al. 2011; Murchison et al. 2012) 
and chromosomal polymorphism (Deakin et al. 2012; Pe- 
arse et al. 2012), DFT cells have been described as a stable, 
clonally evolving cell line. However, our recent studies 
show that this unique cancer is a dynamically evolving obli- 
gate parasite, which uses gene expression alterations (Ujvari 
et al. 2012, 2013), telomere homeostasis (Ujvari et al. 
2012) and epigenetic variations (Ujvari et al. 2013). The 
results from the present study suggest that ploidization 
may offer yet another pathway to which DFTD is able to 
adapt to the ever-changing evolutionary landscape sculp- 
tured by the devils' immune system. Finally, our study is 
the first to show that anthropogenic selection may enhance 
cancer evolution in the wild, and it therefore cautions 
about what measures we employ to try to halt the spread of 
this devastating disease. 
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